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Exciton states in monolayer MoSe2: impact on interband transitions 
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We combine linear and non-linear optical spectroscopy at 4 K with ab initio calculations to study 
the electronic bandstructure of MoSe 2 monolayers. In 1-photon photoluminescence excitation (PLE) 
and reflectivity we measure a separation between the A- and B-exciton emission of 220 meV. In 2- 
photon PLE we detect for the A- and B-exciton the 2p state 180 meV above the respective Is 
state. In second harmonic generation (SHG) spectroscopy we record an enhancement by more than 
2 orders of magnitude of the SHG signal at resonances of the charged exciton and the Is and 2p 
neutral A- and B-exciton. Our post-Density Functional Theory calculations show in the conduction 
band along the K — T direction a local minimum that is energetically and in k-space close to the 
global minimum at the K-point. This has a potentially strong impact on the polarization and energy 
of the excitonic states that govern the interband transitions and marks an important difference to 
M 0 S 2 and WSe 2 monolayers. 


I. INTRODUCTION 


Monolayers (MLs) of the transition metal dichalco- 
genides (TMDCs) M 0 S 2 , MoSe 2 , WS 2 and WSe 2 (abbre¬ 
viated MX 2 ) are semiconductors with a direct bandgap 
in the visible region m- MX 2 MLs are promising ma¬ 
terials for (opto-)electronics [U dj HHH], non-linear op¬ 
tics EHH and for exploring electron k-valley physics 
[IIHIS!. The MX 2 ML materials share common char¬ 
acteristics: (i) Their optical properties are dominated 
by excitons, strongly Goulomb-bound electron hole pairs 
[ISIIMSI]- (ii) In these TMDC MLs crystal inversion 
symmetry breaking together with the strong spin-orbit 
(SO) interaction leads to a coupling of carrier spin and 
k-space valley physics, i.e., the circular polarization 
or a ~) of the absorbed or emitted photon can be directly 
associated with selective carrier excitation in one of the 
two non-equivalent K valleys (K'*' or K“, respectively) 

[laisMszi- 

Using non-resonant laser excitation in ML M 0 S 2 [341- 
[SniEE]) WSe 2 [inilSZ] and WS 2 [33] high values (« 50%) 
for the circular polarization Pc of the stationary photolu¬ 
minescence (PL) corresponding to successful valley polar¬ 
ization have been reported. In contrast, the polarization 
reported for the promising material MoSe 2 [HI UHl SI] 
under similar experimental conditions is surprisingly very 
low (< 5%) [33l444j . Experiments combing PL excitation 
(PLE) and time resolved PL suggest either ultra fast 
polarization relaxation in MoSe 2 in the sub-picosecond 
range or inefficient optical polarization generation due to 
anomalies in the bandstructure to be at the origin of this 
low valley polarization [44] . 

Our target is to investigate how an eventual bandstruc¬ 
ture anomaly (local extremum), can influence the nature 
of the optical transitions. To this aim we combine lin¬ 
ear and non-linear optical spectroscopy at T=4 K with 
calculations beyond standard Density Functional Theory 
(DFT) to study the electronic bandstructure of MoSe 2 
MLs. In 2-photon PLE we detect 180 meV above the 
Is state the 2p state of the A-exciton, well separated in 
energy from the Is B-exciton emission recorded in both 


reflectivity and 1-photon PLE. The energy position of the 
neutral and charged exciton transitions is determined us¬ 
ing second harmonic generation (SHG) spectroscopy, as 
the SHG signal is enhanced by more than 2 orders at 
these particular resonances, important for applications 
in non-linear optics. 

The strong excitonic effects that dominate the 1- and 
2-photon spectroscopy results in the investigated MoSe 2 
MLs have important consequences not just in terms of the 
transition energy [Hj • Due to the small exciton Bohr ra¬ 
dius, the extension in k-space is considerable and the ex¬ 
citon wavefunction will include contributions from states 
far away from the K-point of the Brillouin zone where the 
direct free carrier (electronic) bandgap of the MX 2 MLs 
is situated. The position in energy and in k-space of the 
adjacent local minima in the conduction band (CB) and 
maxima in the valence band (VB), respectively can be 
very different from one MX 2 material to another [JSJ HH] ■ 
Our DFT-GW calculations show in the GB of ML MoSe 2 
a local minimum that is energetically and in k-space close 
to the global minimum at the K-point [47]. We evalu¬ 
ate how the proximity of this local minimum along the 
K — T direction can influence the polarization and en¬ 
ergy of the excitonic states that govern the interband 
transitions, marking an important difference compared 
to M 0 S 2 and WSe 2 monolayers. 

II. SPECTROSCOPY OF EXCITON STATES 

Similar to the hydrogen model, an electron hole pair in 
the TMDC ML interacts through an attractive Coulomb 
potential and will form below the gap a series of excitonic 
Rydberg-like states with definite parity, where ns states 
are of even, and np states are of odd parity (n is an inte¬ 
ger) [3S]. The exciton binding energy E}, of the order of 
0.5 eV in ML MoSe 2 can be estimated by determining the 
free carrier bandgap in scanning tunnelling spectroscopy 
(STS) and subtracting the neutral exciton energy (corre¬ 
sponds to the Is A-exciton state), see Ugeda et al. [25] . 
Here we aim to uncover in addition to the Is exciton 
ground state also higher excited exciton states, which 
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FIG. 1: (a) Conduction (CB) and valence states (VB) in 

a single particle picture, (b) The intensity of the neutral Is 
A-exciton PL (shown as gray solid line) is recorded as a func¬ 
tion of laser energy. In 1-photon PLE (red circles) the Is B- 
exciton state (X||) is identihed, in 2-photon PLE (open black 
squares) the peaks are assigned to the 2p A-exciton (X^) 
and 2p B-exciton (X^). (c) Maxima of the first derivative of 
the reflectivity allow to assign the Is A-exciton and B-exciton 
state energies. 


were shown to dominate the linear and non-linear opti¬ 
cal response in M 0 S 2 and WSe 2 MLs [Ullin]. In Fig.[^ 
we see a typical time-integrated PL spectrum at T=4 K 
for ML MoSe 2 with two prominent peaks (see Appendix 
for information on samples and the experimental set-up). 
The higher energy peak (FWHM=10 meV) at 1.667 eV 
has previously been attributed to the neutral A-exciton 
Xjf [221111111]. At 1.633 eV we record the trion emis¬ 
sion (T) corresponding to a binding energy of 34 meV, 
in agreement with previous measurements [12 IUHlll • 


A. One-Photon PLE experiments 

Here we detect the intensity of the PL emission 
as a function of the laser energy, shown in Fig. [^. We 



FIG. 2: The second harmonic generation (SHG) signal is 
plotted as a function of twice the laser photon energy. The 
local maxima corresponding to the neutral and charged exci- 
ton states are marked (see text). 


find a clear maximum at 1.885 eV, which we attribute to 
the Is B-exciton state. In order to confirm this, we have 
performed reflectivity measurements using a white light 
source. We find in reflectivity two prominent transitions 
in Fig. [TJ:, one at 1.667 eV, which corresponds exactly to 
the neutral A-exciton emission energy measured in PL. 
This indicates negligible localization of excitons in this 
sample (no Stokes shift). The second transition in reflec¬ 
tivity at 1.885 eV gives the position of the Is B-exciton. 
The measured difference between Is A and B-exciton is 
Aab — 220 meV determined from both PL and reflec¬ 
tivity. From our DFT calculations (see section IB for 
details) we obtain a spin splitting of A'^q = 183 meV 
in the valence band. This is in very good agreement 
with the measured value of 180 meV from angle-resolved 
photoemission spectroscopy (ARPES) measurements at 
T=40 K |3]. The CB spin splitting is predicted to have 
the opposite sign compared to the VB, as sketched in 
Fig.[^ [12E2EI]- We obtain from our DFT calculations 
a CB spin splitting of A^q = —30 meV. The measured 
difference of 220 meV is in close agreement with 
the calculated energy difference taking into account the 
VB and CB spin splitting of 183-1-30=213 meV. In addi¬ 
tion, we find in DFT calculations in section |BI| that the 
effective masses in CB and VB are slightly different for 
different spin states, and hence also the exciton binding 
energies of A- and B-excitons are expected to differ and 
contribute to the measured value of Aab- 


B. Two-Photon PLE experiments 

An experimental signature of excited exciton states, 
in analogy with the ns and np states of the hydrogen 
atom, have not been reported yet for MoSe 2 MLs, to the 
best of our knowledge. The excited exciton states and 
the exact exciton binding energy Et in MX 2 MLs are 
currently debated in the literature for WS 2 [12 |12 
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[50] . M 0 S 2 [HI[SI] and WSe 2 nuiMiisni, with differences 
in the reported Et for the same material by factor of 2 , so 
additional experiments are important. In our 2-photon 
PLE measurements we can directly address exciton states 
with p-symmetry and not with s-symmetry as in the 1 - 
photon PLE and reflectivity. In the 2-photon PLE of 
Fig. [ 3 , we find a well defined peak at 1.844 eV, which 
we assign to the 2p A-exciton transition, not obscured 
by the Is B-exciton, which is parity forbidden here. We 
measure an energy difference ls-2p of 180 meV. Using 
a binding energy of the order of E;, « 0.5 eV [H] in a 
simple 2D hydrogen model [15], the energy separation 
ls-2s,2p is expected to be |Et, ~ 440 meV, much larger 
than our measured value of 180 meV. This shows for ML 
MoSe 2 a strong deviation from a simple hydrogenic series, 
observed also in ML WS 2 and WSe 2 , very likely due to 
a strong variation of the effective dielectric constant as a 
function of the spatial extension of the exciton state [52l 
153] . Interestingly we observe a second, clear maximum 
at 2.06 eV, about 175 meV above the B-exciton Is state. 
We tentatively assign this peak in 2-photon PLE to the 
2p B-exciton state. 


C. Second Harmonic Generation Spectroscopy 

In addition to the 2-photon PLE signal we can also plot 
the intensity of the SHG signal as a function of Laser en¬ 
ergy in Fig.[^ The SHG spectroscopy results in this ML 
material without crystal inversion centre also show clear 
resonances at the exciton energies. Interestingly we ob¬ 
serve a strong SHG signal at the transition energy of the 
charged exciton (trion). This hints at a substantial den¬ 
sity of states of this complex. 

The SHG signal at the Is B-exciton is about 120 times 
higher than in the region between 1.7 and 1.8 eV. This 
finding shows that excitonic resonances dominate the 
non-linear optical response of ML MoSe 2 . The SHG sig¬ 
nal is non-zero for all laser energies, the SHG peaks out¬ 
side the exciton resonances between 1.7 and 1.8 eV are 
not clearly visible on the linear scale used in Fig. We 
confirm the energy position of the Is and 2 p exciton lev¬ 
els A and B. Note that resonances of Is exciton states in 
SHG are forbidden if strict electric dipole selection rules 
apply. The prominent Is features in our SHG spectrum 
in Fig. [^ can come from the interplay of electric dipole 
with magnetic dipole transitions HH IMj- In addition 
strict electric dipole selection rules could be slightly re¬ 
laxed if the overall symmetry of the crystal is lowered by 
extrinsic effects such as the substrate. 


III. ELECTRONIC BAND-STRUCTURE 
CALCULATIONS 

In the experiments on ML MoSe 2 we find sharp 
excitonic features that are well defined as for the 
related ML material WSe 2 m- Despite similar exciton 


binding energy and excited state spectrum in ML WSe 2 
and ML MoSe 2 the optical valley polarization using 
non-resonant laser excitation can be generated in the 
former but not in the latter |55l444j . One reason could 
be an anomaly in the bandstructure, which we aim 
to uncover in band structure calculations for MoSe 2 . 
For comparison and to validate our computational 
approach we have performed in parallel calculations 
with the more thoroughly studied ML M 0 S 2 . The 
bandstructure of ML MoSe 2 and M 0 S 2 can be compared 
in Fig. where striking differences in both valence 
and conduction band appear. We extract the effective 
carrier masses, the band gaps and the exciton binding 
energies. We discuss the most important features for 
optical transitions in detail, in particular the competition 
between direct and indirect exciton states in MoSe 2 MLs. 


A. Quasi-Particle Band Structure 

The comparison of quasi-particle band structure cal¬ 
culated by DFT (see Appendix for computational de¬ 
tails) of ML M 0 S 2 and MoSe 2 provides interesting in¬ 
sights: when replacing S with Se, the spin-splitting is en¬ 
hanced. For M 0 S 2 the valence and conduction state the 
spin-splittings at the K+ point are = 143 meV and 
AgQ = —13 meV, respectively. For ML MoSe 2 we find 
AgQ = 183 meV and A^q = —30 meV. The direct band 
gap values (neglecting excitonic effects) at K+, depicted 
by arrows in the left panels of Fig. are 2.13 and 2.31 
eV for MoSe 2 and M 0 S 2 , respectively. The latter value is 
lower than previous theoretical studies as energies from 
2.41 to 2.97 eV have been reported [HI IHl UHl i55H55] . 
It is important to take into account that the uncertainty 
in this type of band gap calculations is in the hundreds 
of meV range, depending on the computational settings. 
The obtained results critically depend on the choice of 
the GW methods used (self-consistent or partially self- 
consistent scheme), the number of unoccupied states in¬ 
cluded, the vacuum height and the k-point sampling. 
Interestingly using exactly the same computational set¬ 
tings, our MoSe 2 direct band gap estimate is closer to 
recent reports: 2.33 eV for a Go Wo calculation in Refi59l 
when Ugeda et al [H] propose 2.26 eV for a free standing 
MoSe 2 ML. On the experimental side, using STS tech¬ 
niques the band gap of MoSe 2 ML is 2.18 ± 0.04 eV [37] 
but it includes substrate screening effects. 

If one compares the full band-structure of the two 
material systems another striking difference appears: in 
M 0 S 2 ML, the topmost valence bands in F and K+ are 
only separated by 65 meV. This feature appears to be 
essential in the mechanism of direct to indirect band gap 
transition and its evolution with layer thickness (IML, 
2ML ...) [60] and is expected to have a strong impact on 
the optically generated valley polarization as a function 
of laser energy [3^ [Ml . The same trend for the gap evo¬ 
lution with layer thickness applies to MoSe 2 [3], albeit 
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FIG. 3: GoVFo band structure for freestanding MoSe 2 ML (top left). Zoom at K point and along K-F direction (top right). 
M 0 S 2 ML results are in the bottom panels. Zero energy represents the Valence Band Maximum. The green shaded area is a 
schematic representation of the extension in k-space of the Is A-exciton wavefunction, see. Eq. (4) in the text. 


with a much larger T to K"*" energy separation of 485 
meV for IML. 

For ML MoSe 2 our calculations show a remarkable 
anomaly in the CB, which could play a key role for the 
valley polarization dynamics: In ML MoSe 2 there is a 
small energy difference of only 46 meV between the CB 
minimum in K+ and in A, as can be clearly seen in Fig.[^ 
The corresponding value for M 0 S 2 is 345 meV and there¬ 
fore not in competition with direct optical transitions at 
the K points. We note another important difference: The 
A position along the K—F line for M 0 S 2 lies exactly at 
the (1/6,1/6, 0) coordinates. In contrast, for MoSe 2 it 
is slightly shifted towards the K point. Interestingly the 
spin-up and spin-down minima are not positioned at the 
same k-value. Note that the appearance of this local min¬ 
imum has also been predicted recently in m- 

In the following we discuss how the proximity of the 
CB minimum at A can potentially influence the optical 
transitions. Excitons with large binding energies in ML 
MoSe 2 have a very small Bohr radius and as a result ex¬ 
tend in k-space well beyond the K-point (see Fig. |^. 
To allow for a more quantitative discussion, we need 
to extract the effective carrier masses from our calcu¬ 
lated band structure. At VB and CB extrema, effective 


masses for holes and electrons can be extracted with a 
simple quadratic fitting procedure. The extracted mass 
values will depend on the selected k-value range around 
the extrema. For consistency, we have used the inter¬ 
val [—fco/2, +ko/2], with ko being defined as 1/(2 x o^^), 
where is the Bohr radius of the exciton ground state 
in a Wannier-Mott picture of electron-hole pair [27]. The 
computed effective mass values are given in Table [Tj 
For M 0 S 2 , the rUc and values are only slightly dif¬ 
ferent from previous studies |20l [55] [58], the difference 
being reduced in the electron-hole pair effective mass 

p. = - -—— values. There is a certain spread in values, 

nic + rriv 

especially for the M 0 S 2 hole effective mass with the value 
provided in Ref. |57| or for MoSe 2 values [SS] . We obtain 
smaller values, but we attribute these differences to the 
choice of the extension in k-space around the extrema. 
Indeed, if one takes a fitting interval of length 2/co the ef¬ 
fective masses are significantly increased by 30% in both 
material systems. Interestingly spin-up and spin-down 
excitons have different effective masses at K+: 0.21 mg 
vs 0.24 mo for and respectively for M 0 S 2 , 0.25 mg 
and 0.28 mg in MoSe 2 , due to larger values of the spin- 
down rric and rriy, in both systems. This confirms that 
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the exciton binding energies for A- and B-excitons should 
be slightly different, which will in turn influence the en¬ 
ergy difference between the corresponding optical 

transitions measured in Fig. 

On top of Go Wo calculations, exciton binding energies 
Ah, given in Table |lj have been extracted from the imag¬ 
inary part of transverse dielectric constant after solving 
the Bethe-Salpeter Equation (BSE). Our 0.58 eV esti¬ 
mate is in reasonable agreement with previous theoreti¬ 
cal studies for the M 0 S 2 A-exciton binding energy: values 
from 0.55 [52], around 0.6 [IHlIiHI, or around 0.9 ^51 [57] 
and up to 1.1 eVmilSSj have been reported. Please note 
the computational details are extremely different for each 
of these values. It is therefore difficult to extract precise 
trends. Alternatively, the spread in values should rather 
be used as a reasonable error bar for this type of calcu¬ 
lations. For MoSe 2 theoretical values for Af, of 0.9 [55] . 
0.78 [21] and 0.65 eV [25] can be found in the recent liter¬ 
ature. This indicates that the calculated exciton binding 
energy in ML MoSe 2 is smaller than in M 0 S 2 , as observed 
experimentally. 


TABLE I: Free carrier band-gap Eg, exciton binding energy 
El,, CB effective masses ra^ and VB effective masses m„ (here 
given as positive values as in the hole picture) and the cor¬ 
responding relative static dielectric constant and exciton 
Bohr radius for M 0 S 2 and MoSe 2 Monolayer (see text). 



E, 

El, 

rric 

rriv 

tr 

O-B 

in 

eV 

eV 

rrie 

me 

- 

A 

M 0 S 2 

2.31 

0.58 

0.40 

0.46 

4.47 

2.77 

MoSe 2 

2.13 

0.51 

0.49 

0.52 

5.17 

2.40 


B. Exciton States 


Our target is now to estimate how far the exciton state 
extends in k-space around the K+ point. This allows 
us to estimate the possible impact of the proximity of 
the local CB minimum at A on the optical transitions 
and the comparatively weak valley polarization in ML 
MoSe 2 . This order of magnitude discussion is graphically 
represented in the right column of Fig. [^ Starting from 
the standard 2D Wannier-Mott model [48] , the exciton 
ground state energy writes simply as: 


Eb 




( 1 ) 


ql 


with e*^ = — = ———. An estimate of the correspond¬ 
ed dTregCd 

ing exciton Bohr radius is thus given by 




( 2 ) 


Combining these equations allows for the calculation of 
the relative dielectric constant, knowing the binding en¬ 
ergy Eb and the reduced mass of the electron-hole pair, 
extracted from GoWq-I-BSE calculations. It yields values 
of 4.47 and 5.17 for M 0 S 2 and MoSe 2 , respectively. The 
main differences from previous theoretical results arise 
from our choice of a 2D model, the use of our own com¬ 
puted Af, and fx. These values are usually larger than the 
estimates found in Ref. [2ni ESI using K = y'ej_e|| (3.44) or 
by our own direct calculation of the static dielectric con¬ 
stant at the DFT level including local field effect (2.61) 
for M 0 S 2 . Using the same k approximation we find 2.80 
for MoSe 2 in comparison with 5.17 of Table [l] 

Continuing our discussion based on a simple hydrogen¬ 
like model, the exciton ground-state wave function can be 
written as: 


01s w 



(3) 


where Og = and r is the coordinate of the relative 
electron-hole motion. The corresponding Fourier trans¬ 
form of this ground state provides a rough estimation of 
how far around the K+ point the exciton state is spread 
in k-space. 


01s(9) 


(l + (2 


q 


2D\ 


3/2 ’ 


(4) 


Our representation of 4>is{q) for ML MoSe 2 in Fig. [^sug¬ 
gest the following scenario: The close proximity of the A 
conduction band minimum could be at the origin of the 
low valley polarization degree of the A-exciton Is state 
luminescence, as contributions to the exciton wavefunc- 
tions away from the points do not obey the strict 
chiral valley selections rules. It is important to underline 
in Fig.[^the striking difference with ML M 0 S 2 , where the 
exciton state is mainly build with K'*" electronic states as 
the contribution from the states around the A point at 
much higher energy is negligible. 


IV. CONCLUSIONS 

We probe the exciton states in ML MoSe 2 in 1- and 2- 
photon PLE. We see a clear signature of the 2p state of 
the A- and B-exciton about 180 meV above the respec¬ 
tive Is exciton state in 2-photon PLE. Our post-DFT 
calculations reveal in the MoSe 2 ML conduction band a 
local minimum at the A point only 46 meV above the 
global minimum at the K-points. As excitons with large 
binding energies (Af, « 0.5 eV) involve k-states far away 
from the K-point, the contribution of states around the 
A point to optical transitions is possible. As here the 
chiral optical selections rules are not applicable, these 
states could contribute to the overall low valley polar¬ 
ization reported for MoSe 2 To further test this 

hypothesis, experiments that modify the band structure 
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via strain tuning [631164] would provide useful informa¬ 
tion. Also a systematic study of valley polarization in 
MoS 2 Se 2 alloy monolayers [SS] will help clarifying, why 
M 0 S 2 shows strong and MoSe 2 shows very weak valley 
polarization in non-resonant PL experiments. 
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V. APPENDIX 

A. Samples and Experimental Set-up 

MoSe 2 ML flakes are obtained by micro-mechanical 
cleavage of a bulk MoSe 2 crystal on an Si 02 /Si sub¬ 
strates using viscoelastic stamping [B5]. The ML region 
is identified by optical contrast and very clearly in PL 
spectroscopy. Experiments at T=4 K are carried out in 
a confocal microscope optimized for polarized PL exper¬ 
iments m- The MoSe 2 ML is excited by picosecond 
pulses generated by a tunable frequency-doubled optical 
parametric oscillator (OPO) synchronously pumped by a 
mode-locked Ti:Sa laser. The typical pulse and spectral 
width are 1.6 ps and 3 meV respectively; the repetition 
rate is 80 MHz. The laser average power is tunable from 
2 to 200 /rW. The detection spot diameter is « l^m, i.e. 
considerably smaller than the ML diameter. For time in¬ 
tegrated experiments, the PL emission is dispersed in a 
spectrometer and detected with a Si-GGD camera. 


B. Computational Details 


The atomic structures, the quasi-particle band struc¬ 
tures and optical spectra are obtained from DFT calcu¬ 
lations using the VASP package [Ml IM] ■ PBE functional 
[70] is used as approximation of the exchange-correlation 
term. It uses the plane-augmented wave scheme 
to treat core electrons. Fourteen and six electrons for 
Mo and S, Se respectively are explicitly included in the 
valence. All atoms are allowed to relax with a force 
convergence criterion below 0.005 eV/A. After primi¬ 
tive cell relaxation, the optimized lattice parameters are 
3.22 and 3.32 A for M 0 S 2 and MoSe 2 respectively, these 
values being in good agreement ( 1%) with previous 
studies [Ml ES] and slightly larger than the bulk exper¬ 
imental values. To sample the Brillouin zone a grid of 
12 x12x1 k-points has been used, in conjunction with a 
vacuum height of 17 A, to take benefit of error’s can¬ 
cellation in the band gap estimates [M]- This provides 
exciton binding energies in reasonable agreement with 
experiments as suggested in different works [28l [56] . A 
gaussian smearing with a width of 0.05 eV is used for 
partial occupancies, when a tight electronic minimization 
tolerance of 10“® eV is set to determine with a good pre¬ 
cision the corresponding derivative of the orbitals with 
respect to k needed in quasi-particle band structure cal¬ 
culations. Spin-orbit coupling was also included non-self- 
consistently to determine eigenvalues and wave functions 
as input for the full-frequency-dependent GW calcula¬ 
tions m performed at the GqWq level. The total number 
of states included in the GW procedure is set to 600, after 
a careful check of the direct band gap convergence smaller 
than 0.1 eV. We have used the WANNIER90 code [74] 
and the VASP2WANNIER90 interface [7S] to interpolate 
the band structures on a finer grid. Optical absorption 
spectra have been calculated using Bethe-Salpeter Equa¬ 
tion in the Tamm-Dancoff approximation, including the 
six highest valence bands and the eight lowest conduction 
bands [T5] . 
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